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Abstract. Recently, Huang and Lin suggested a combination of two successful! 
mean-field theories, the 2-cluster approximation and paradisical mean-field, for 
the Nagel-Schreckenberg cellular automaton model of traffic fiow. We argue that 
this new approximation is inconsistent since it violates the Kolmogorov conditions. 



In a recent work Huang and Lin have studied a combination of two mean-field 
theories, the 2-cluster approximation ^, Q and paradisical mean- field (PMF) |,|, 
for the Nagel-Schreckenberg cellular automaton model of traffic flow |^ (for a review, 
see 0). The suggested combined theory yields results for the flow-density relation 
which are worse compared to Monte Carlo simulations than those of each individual 
approximation alone. The authors concluded that the success of paradisical mean- 
field theory is accidental and cannot be improved systematically, in contrast to the 
cluster approximation. In this comment we will show that the combined theory as 
suggested by Huang and Lin is inconsistent since it violates the elementary Kolmogorov 
conditions. This also invalidates the conclusions about the success of paradisical mean- 
field theory. 

In the 2-cluster approximation the probability P(ti, . . . , r^) to find the system in 
the state (ti, . . . ,tl) is factorized into 2-site probabilities Pt,t''- 

P{t\ , . . . , T/^) OC Pri ,T2 ,T3 ' ' ' -^Tl-I.Tl ,ti ■ (1) 

Here we assume translational invariance so that the probabilities Pr.r' do not depend 
on the position. 

An important consequence are the so-called Kolmogorov consistency conditions 
which for the 2-cluster approximation read 

, Pri,T2 = Pti = ^ ^ Pt2.ti (2) 



r2 T2 



where Pr is the probability to find a cell in state r. Especially for t — denoting an 
empty cell in the notation of , we have Px = 1 — p where p is the total density of 
cars. Since also Pi + P2 = p, instead of eq. (2) in [Q one has more precisely 

P = Pix + Pii + P12 + P2x + P21 + P22 (3) 
= Pxi + Pu + P21 + Px2 + P12 + P22 (4) 
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and 

l-p=Pxx+ Pxl + Px2 (5) 
= Pxx+Plx+P2x- (6) 

We have analyzed numerically the equations given in the Appendix of 0| and 
found that these consistency conditions are violated. In fact we found that the 
results depend strongly on the iteration scheme used. If one adjusts the normalization 
M during the iteration to ensure the normalization of probabilities (eq. (1) of 
then the numerical solution converges to p = 0, even for initial values with p ^ Q. 
On the other hand, if one uses eq. (1) of 0] to eliminate one of the equations in 
the Appendix, the normalization condition is always satisfied during iteration. One 
then can use M to fix the density p according to eq. (2) of Q. In this case we 
find that the Kolmogorov consistency conditions (||)-(|^) are violated and in general 
^^(PiT- + P2t) 7^ P 7^ TliTiPri + Pt2)- It even seems that no solution exists, at least 
in certain density regimes. 

In order to exclude the possibility that these problems are related to a mistake 
or typo in the set of equations given in we have done an independent calculation. 
However, we find it much more convenient not to use the update ordering R2-R3- 
R4-R1 that was used in ||, ^ and adopted in [0. Instead the original update order 
R1-R2-R3-R4 |^ is preferable since it allows a much easier identification of the GoE 
states and avoids the introduction of the weighting W. It also turns out that this 
approach is numerically much more stable since no terms proportional to (where 
M is the normalization) appear. For our system of equations we find the same type 
of behaviour as described above for the equations of Q . 

The reason for this failure is that the modified equations after elimination of 
the GoE states do not automatically guarantee the conservation of cars, in contrast 
to the original 2-cluster equations. It is not possible to satisfy the normalization 
condition r2^'ri.^2 ^ 1 the Kolmogorov conditions (j^) by just introducing 
one normalization factor M . One could try to use normalizitions A/^.t' depending on 
the state, but this would make the system of equations even more complicated than 
the 3-cluster approximation ||, ^ which already treats most GoE states correctly. 

Finally we want to comment on the conclusions of We believe that the 
combined approach - if successful - must be considered as an improvement of the 
2-cluster approach, not of PMF as suggested in [0. If it would be possible to carry out 
such an approximation consistently it would most certainly lead to an improvement 
of the 2-cluster results. In fact the main point of the PMF is not the quantitative 
agreement with simulations but the fact that it allows to determine the origin of 
correlations and thus to improve the understanding of the underlying physics. This 
was already emphasized in 1 1, Q . 
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